Photoemission spectrum and effect of inhomogeneous pairing fluctuations in the 
BCS-BEC crossover regime of an ultracold Fermi gas 
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We investigate the pliotoemission-type spectrum in a cold Fermi gas which was recently measured 
by JILA group [J. T. Stewart et ai, Nature 454, 744 (2008)]. This quantity gives us very useful 
information about single-particle properties in the BCS-BEC crossover. In this letter, including 
pairing fluctuations within a T-matrix theory, as well as effects of a harmonic trap within the local 
density approximation, we show that spatially inhomogeneous pairing fluctuations due to the trap 
potential is an important key to understand the observed spectrum. In the crossover region, while 
strong pairing fluctuations lead to the so-called pseudogap phenomenon in the trap center, such 
strong-coupling effects are found to be weak around the edge of the gas. Our results including this 
effect are shown to agree well with the recent photoemission data by JILA group. 

PACS numbers: 03.75.Ss,05.30.Fk,67.85.-d 



The recent photoemission-type experiment developed 
by JILA group [ij provides a powerful method to study 
microscopic properties of cold Fermi gases. This ex- 
periment is an analogue of angle-resolved photoemission 
spectroscopy (ARPES)i2], which has been extensively 
applied in condensed matter physics. Using this tech- 
nique, one can probe single-particle excitations that al- 
low us to investigate many-body effects in the BCS-BEC 
crossover .3- 14]. Indeed, the observed spectra exhibit 
dramatic change in the crossover region [l|. 

One typical many-body effect on single-particle excita- 
tions expected in the BCS-BEC crossover regime of a cold 
Fermi gas is the pseudogap effect. In this phenomenon, 
preformed pairs cause a gap-like structure in the den- 
sity of states (DOS) even above the superfluid phase 
transition temperature Tc- The pseudogap has been ob- 
served in the underdoped regime of high-Tc cupratesfTsj. 
However, as the origin of this phenomenon in high-Tc 
cuprates, in addition to pairing fluctuations, various pos- 
sibilities have been proposed, such as spin fluctuations 
and a hidden order. Since the BCS-BEC crossover in a 
cold Fermi gas is dominated by pairing fluctuations, the 
study of pseudogap in this system would be also useful 
for clariiying the validity of the pseudogap mechanism 
based on preformed pairs in high-Tc cuprates. Since the 
pseudogap appears in single-particle excitations, the pho- 
toemission experiment would be useful for this purpose. 

In considering the pseudogap effect in cold Fermi gases, 
one should note that the presence of a trap naturally 
leads to spatially inhomogeneous pairing fluctuations. 
While the pseudogap structure in DOS is expected to 
be remarkable in the trap center, such a many-body ef- 
fect may be weak around the edge of gas cloud. Since the 
pseudogap in high-Tc cuprates is a uniform phenomenon, 
this inhomogeneous pseudogap effect is unique to trapped 
Fermi gases. Indeed, the double-peak structure in the 
photoemission spectrum observed in the strong-coupling 
BEC regime of Fermi gas[l| seems difficult to explain 
as far as a simple uniform system is considered (l6|. 



In this paper, we study pseudogap phenomena in a 
trapped Fermi gas above Tc, addressing the recent pho- 
toemission experiment by JILA group [l|. Extending our 
previous paper for a uniform Fermi gas to include effects 
of a harmonic trap within the local density approxima- 
tion (LDA), we calculate the local DOS, as well as the 
local spectral weight (SW), over the entire BCS-BEC 
crossover region within the T-matrix approximation in 
terms of pairing fluctuations [l6l - l20| . We clarify how the 
inhomogeneous pseudogap phenomenon appears in these 
quantities. 

Recently, Refs. [2l], have studied the photoemis- 
sion spectra using phenomenological theories with the 
BCS ansatz. While they capture some features of experi- 
mental results in the unitarity limit, they cannot explain 
the observed spectrum in the BEC regime, which con- 
sists of upper sharp and lower broad peaks. In contrast, 
including the inhomogeneous strong-coupling effect, we 
show that our ab initio calculation of the photoemission 
spectra can naturally explain the experiments 1] from the 
weak-coupling BCS to the strong-coupling BEC regimes, 
without introducing any free parameters. 

We consider a two-component Fermi gas in a harmonic 
trap. Assuming a broad Feshbach resonance, we employ 
the ordinary BCS model, described by the Hamiltonian, 



^ E '^p+q/2t'^-p+q/2iV+q/2iV+q/2t- 
<1,P,P' 



(1) 



Here, 



is the annihilation operator of a Fermi atom 



with pseudospin a =t, i, and the kinetic energy = 
£p — /i is measured from the chemical potential /i, where 
£p = /2m and m is the atomic mass. The pairing inter- 
action —Ui< 0) is assumed to be tunable by a Feshbach 
resonance |23|. which is related to the s-wave scattering 
length as as^ inas/m ^ -U/[l ~ UJ2pl/{2ep)]. 
Within a LDA, effects of a trap are conveniently in- 
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FIG. 1: Calculated Tc as a function of the inverse scattering 
. The inset shows fj, at Tc- kp is the Fermi momentum, 
and £f is the Fermi energy. We use these Tc and p in FiES. I2l5l 



corporated into the theory by replacing the chemical 
potential ^ by fi{r) = — F(r)[2^, where V{r) — 
is a harmonic potential (where Wtr is a trap 
frequency). The LDA single-particle Green's function 
is given by Gpiiu^.r) = l/[iw„ - ^p(r) - Ep(zw„,r)], 
where ^p(r) = Ep — mC'')! ^^d aj„ is the fermion Mat- 
subara frequency. The LDA self-energy I]p(iaj„,r) in- 
volves effects of pairing fluctuations within the T-matrix 
approximation [16|-|20|, which is given by I]p(iaj„,r) — 
'^I]q,i,.„ rq(w„,r)G°_p(w„ - iun.r). Here, m„ is 
the boson Matsubara frequency, and Gp(iaj„,r) — 
— £,p{t)\ is the LDA free fermion propagator. 



-U/[l — Uliq{ivn.,r)] is the particle- 



particle scattering matrix within the T-matrix approx- 
imation, where Xlq[wn,r) = ^ Sp.ii.^^ ^p+ q/2( ^^" + 
«cj„, r)G'^p^q^2(~*'^ni ^) is the pair propagator [lOj. 

In this paper, we focus on the case at Tc to compare 
with the experiments in Ref. The region above Tc 
will be discussed elsewhere. In LDA, the superfluid phase 
transition is determined by the Thouless criterion at the 
trap centerd^, i.e., Tq=o{iVn = 0, r = 0)"^ = o[Tot. 
Strong-coupling effects on /i arc included by solving this 
Tc-equation together with the equation for the total num- 
ber of Fermi atoms, N = 2T J dr Ep,»c.„ Gp(«w„, r)e"^"^ . 

Figure [T] shows the calculated Tc and /i as functions of 
{kpas)~^, where fci? is the Fermi momentum defined from 
the Fermi energy ei? — (3A^)^/^a;tr = kp/{2m). In Fig.[Tl 
Tc in a trapped gas is found to be higher than that in a 
uniform gas in terms of e_F[2l,[23]- Using them, we calcu- 
late the local SW A{p,uj,r) and local DOS p{io,r) from 
the analytic continued Green's function as, respectively. 



A{p, 



— Im[Gp(ia;„ 

TT 



+ iS,r)l (2) 



(3) 



The photoemission spectrumQ can be calculated in 
the same way as the rf-tunneling current spectroscopy [26, 
m, 113 , where atoms in one of the two pseudospin states 
(= I t)) are outcoupled into an unoccupied state |3) 
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FIG. 2: (color onhne). (a)-(c) Intensity of local SW yl(p, r) 
(units of Ep^). (d) local DOS p{u;,r). We set T = Tc and 
(fcpas)"^ = 0. The dotted line in (a) is the BCS-like quasi- 
particle spectrum with Apg = 1.06. 



I ti 4-)) by an applied rf pulse. Since the final-state in- 
teraction can be safely neglected in ""^K Fermi gas 1, 21], 
the Hamiltonian for |3) is simply given by H3 = X]p[^p + 
^3 ~ M3('")]^p^p- Here, bp describes the third pseudospin 
state 1 3), and W3 is the energy difference between | t) 
and 1 3). The chemical potential iJ-sir) = /is — V{r) in- 
volves trap effects within LDA. The transition from _H 
to |3) is induced by the tunneling Hamiltonian[26, 211, 23| 
Ht = X]fc[s *"^*^fc+qj^'^fet + H.C.], where tp is a trans- 
fer matrix element, and and ujl are the momen- 
tum and frequency of the rf-pulse, respectively. The 
photoemission spectrum is obtained from rf-tunneling 
current I{il,r) from | f) to |3), where = col — UI3 
is the rf-detuning. Within the linear response theory, 
we obtain I{n,r) - +^ +'M\.5t' 

where J(i) = 

H.c.Je"'*^''^"'"^^'* is the tunneling current operator 
in the Heisenberg representation, and Hxit) = 
gi[H+H3)t^^^-r{H+H3)t _ -y^g ^-^^^ obtain the rf-tunnefing 

current /(f2,r) — '^pI(p,^Tf')^ where the momentum- 
resolved photoemission spectrum /(p, Q,, r) has the form 



. -^^_^dt' {[J{r,t),HT(r,t')])e'' 



I{p,n,r) = 2nt%Aip,^pir) ~ n,r)fiCpir) - n). (4) 

Here, f{il) is the Fermi distribution function. In Eq. (U), 
we have assumed that the momentum of rf-photon is 
neghgible and |3) is initially empty (/(Cp — f^sir)) = 0). 

In a uniform Fermi gas, the photoemission spectrum 
is related to SW of occupied states as I{p, — ^p — 
uj) — 2nt'j^A{p,uj)f{uj). In particular, it is equal to SW 
below a; = at T = 0. When T > 0, thermally excited 
quasiparticles also contribute to the spectrum, so that 
I{p, — > ^p — w) becomes finite even when cj > 0. 

In the photoemission experiment [l[, since the rf-pulse 
is applied to the whole gas cloud, the observed spec- 
trum involves contributions from all spatial regions of the 
cloud. To include this situation, we should take spatial 
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FIG. 3: (color online). Intensity of local SW A(p,Lo,r) (units 
of e;i). (a)-(c) (fcpa^)-! = -1. (d)-(f) (fepo.)-^ = +1. We 
set T = Tc. 

average of Eq. dU as 

(5) 

Here, = AirRp/S, where Rp = \/2^/ (muj^j.) is the 
Thomas- Fermi radius [33|. We emphasize that this gives 
a proper definition for spatially averaged photoemission 
spectrum. For later convenience, we define the averaged 
occupied SW and DOS by, respectively, 

A{p,Uj)f{Lj) = 4vo(p, n^^p^ Lj)/{27Ttl), (6) 

piLu)f{Lj) = J2 ^avc(p, n^^p- L0)/{27Tt%). (7) 
P 

To see the basic characters of Eqs. dH) and (O, it is 
helpful to consider a free Fermi gas at T = 0. In this 
case, Eqs. (|6]) and ^ reduce to, respectively, 

A{p,co)f{oj) = - ^p)e{-oj), (8) 

(9) 

In the former, the peak position gives the one-particle en- 
ergy ^p. In the latter, DOS in a uniform gas (cx y/u + fi) 
is modified by the factor jw^/^. 

We now show our numerical results. Figure[2]shows the 
local SW A{p, uj, r) and DOS p{uj, r) at Tc in the unitar- 
ity limit {{kpas)~^ = 0). In the trap center (panel (a)), 
a clear pseudogap structure exists, i.e., two prominent 
peaks appear along the particle branch and hole branch of 

the BCS-like quasiparticle spectrum oj = ±./^p-|-Apg, 
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FIG. 4: (color online). Intensity of averaged occupied SW 
A{p,uj)f{uj) (units of l/2m) at Tc. The values of pairing 
interaction {kpCLs)'^ are (a) -1, (b) 0, (c) 0.4, and (d) 1. The 
dashed line indicates the spectrum of a free Fermi gas ^p. The 
solid line is the lower peak position of intensity. Symbols are 
corresponding experimental data [ij. 



where the superfluid gap A is replaced by the pseudo- 
gap Apg. The origin of Ap,^ is a particle- hole coupling 
by pairing fluctuations fl6l-l2q| . Since pairing fluctuations 
also induce a finite lifetime of quasiparticle excitations, 
the particle and hole branches in Fig. [UJa) have finite 
widths, which is in contrast to the mean-field BCS case, 
where both branches appear as (5-functional peaks. We 
note that the deviation of the lower peak from the BCS- 
like quasiparticle spectrum is considered due to the pres- 
ence of excited pairs [H, [13] . 

Pairing fluctuations are weak around the edge of the 
gas due to low particle density. Thus, the pseudogap in 
the local SW gradually disappears, as one leaves from 
the trap center. (See Figs. Hfb) and (c).) In panel (c), 
a single sharp peak line only exists near the one-particle 
energy of a free Fermi gas ^p(i?F)- Namely, the gas is 
pseudogapped in the center of the trap, while it is nearly 
non-interacting on its edges. 

These inhomogeneous features can be also seen in lo- 
cal DOS. As shown in Fig.[2ljd), while a large dip struc- 
ture (which is a characteristic pseudogap effect in DOS) 
appears around cj = in the trap center, local DOS 
is almost equal to that of a free Fermi gas (p(a;,r) — 
(TO^/2/(\/27r2))^w -t- p,{r)) when r = Rp. Since spatial 
inhomogeneity is unique to trapped Fermi gases, the ob- 
servation of local SW and DOS by using the tomographic 
techniques [Sll. [3^ would be interesting. 

We note that the inhomogeneous pseudogap structure 
depends on the strength of pairing interaction. When 
pairing fluctuations are weak in the BCS regime, the 
double-peak structure in A(p, aj,r) soon disappears, as 
one leaves from the trap center. (See Figs. [3Ua)-(c).) In 
contrast, Figs.[3][d)-(f) show that the pseudogap features 
persist even near the edge of the gas in the BEC regime. 

Figures |4] and [SJ respectively, show the averaged oc- 
cupied SW and DOS at T^. In the weak-coupling BCS 
regime, these quantities are expected to be close to the 
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FIG. 5: (color online). Averaged occupied DOS p{uu)f{ijj) 
at Tc. The inset shows a magnification of the curve for 
{kpas)~^ = 1. Circles and crosses show the experimental 
data for (kpas)'^ = and 1, respectively 0]. The data are 
fitted so that the peak value coincides with that of theory 
curves. The data agree well with the theory curves. 

free Fermi gas results in Eqs. dS]) and (0). Indeed, the 
position of the peak line in Fig. lUJa) is almost at the 
single-particle energy of a free fermion uj = (< 0), 
and the averaged occupied DOS in Fig. [5] is close to 
DOS for a free Fermi gas multiplied by when 
(fepOs)^^ — ~1- However, we also find that the peak line 
in Fig. m^a) is slightly below the curve of = ^p, which 
is a signature of pseudogap effect near the trap center. 
Namely, the pseudogap Apg lowers the hole branch as 

~l^p(^)l ~^ — ^i^p(r) + A^g, and this effect still remains 
even after spaliial average. 

In the crossover region, pairing fluctuations are strong 
around the trap center. In the unitarity limit, while the 
region around the e dge of the ga s still gives a sharp peak 
line at a; ~ in A{p,uj)f{uj), the pseudogap in SW 
around r — causes the broadening of the lower part 
of the peak line around p/kp ~ 0.5 shown in Fig. Hl^b). 
In addition, short lifetime of quasiparticle excitations by 
strong pairing fluctuations around the trap center also 
causes the broadening of p{uj)f{Lu), as shown in Fig. \E\ 



We also note that one can see the back-bending of the 
lower peak position in Figs. 4(b)-(d) originating from the 
lower branch shown in Fig. 2(a), which well agrees with 
the experimental dataflj. 

In the BEG regime, the pseudogap features persist to 
the edge of the gas. Thus, the double-peak structure in 
the local SW (See Figs. [3Ud)^(f).) is not smeared out by 
spatial average, as shown in Figs. |4ljc)^(d). In particu- 
lar, the double-peak structure in panel (d), consisting of 
upper sharp peak and lower broad peak, agrees well with 
the photoemission experiment by JILA group[lJ. We em- 
phasize that this spectral structure cannot be explained 
within the previous phenomenological theories [U [2^ . 
Such a double-peak structure also appears in Fig. [5l 
which is also consistent with the experiment 

In the BEG limit, the spectral weight only has the up- 
per sharp peak describing the dissociation of two-body 
bound molecules. On the other hand, as discussed in 
Ref. [l^ , the lower broad peak in SW in the BEG regime 
is an evidence of many-body character of paired atoms. 
Thus, the double-peak structure in Figs. |4l[c)~(d) may 
be understood as a result of the fact that, in the BGS- 
BEG crossover, the character of fermion pair continu- 
ously changes from the many-body bound state associ- 
ated with the Gooper instability to the two-body bound 
state where the Fermi surface is not necessary. 

To summarize, we have discussed inhomogeneous pseu- 
dogap effect in the BGS-BEG crossover regime of a 
trapped Fermi gas. Including this effect, we showed 
that calculated spatially averaged photoemission spec- 
trum agrees well with the recent experiment by JILA 
group Since spatial inhomogeneity by a trap is in- 
evitable in a cold Fermi gas, our results would be useful 
for clarifying strong-coupling phenomena of this system 
over the entire BGS-BEG crossover region. 

We wish to thank A. Griffin for useful discussions. We 
acknowledge J. P. Gaebler and D. S. Jin for providing us 
with their experimental data. 
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